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OO \ The increasing cosmic ray statistics collected by present experiments and 
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the future prospects with new large arrays demand accurate calculations of the 



(N 
> 

ff^ , extensive air shower parameters (EAS). The energy of the primary particle is 

(N 

\ estimated by ground arrays fitting a lateral distribution function (LDF) to the 

' particle densities at a given observing level. However, the lack of appropriate 



parameterization for these distributions, able to reproduce the data collected 



O \ from all arrival directions, makes it difficult the experimental analysis. We 

+^ ■ 

^ \ propose a method to parametrize particle density distributions of EAS at any 



incident zenith angle. Starting from analytical LDF for vertical showers we 



H ' present a detailed study of the atmospheric depth dependence of the shower 

parameters. The results obtained are used to calculate the corresponding 
LDF for non-vertical showers including for the first time both, geometrical 
and atmospheric attenuation effects. We check the method analysing electron 
and muon LDF generated by Monte Carlo simulations from incident cosmic 
ray particles at different zenith angles. A comparison of the proposed LDF 
with experimental results as well as MC data including detector effects is also 
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presented. 
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I. INTRODUCTION 



The observation of cosmic ray particles with energies well beyond the predicted Greisen- 
Zatsepin-Kuzmin cutoff enhanced by the discrepancies in the cosmic ray spectrum 

recently reported by HiRes and AGASA experiment is pushing astroparticle physics 
into a period of fast development both theoretically and experimentally. 

The analysis of the extensive air showers (EAS) produced when the high energy particle 
enters in the atmosphere and interacts with the air molecules is crucial to know the charac- 
teristics of the ultra high energy cosmic rays (UHECR). On the other hand, studies of the 
high energy interaction taking place in the atmosphere could yield unique information for 
elementary particle physics. 

The interest in EAS studies has been revived mainly because improvements achieved 
with new large arrays and the increasing statistics collected demand accurate calculations 
of the shower parameters. A very important shower observable is the spatial distribution 
of particles at a given observing level because ground array techniques use the particle 
density distributions of EAS to estimate the energy of the primary particle, fitting a lateral 
distribution function (LDF) to the observed particle densities 0] . There exist several LDF 
parameterizations well suited to describe density distributions for primary incident angles 
smaller than 45 degrees, which are currently used by cosmic ray experiments. 

High zenith angle showers induced by protons, heavy nuclei or photons have been sub- 
ject of recent studies as an important tool to obtain relevant information about CR mass 
composition [^H. Very inclined showers are also of great importance since they are the 
main background for neutrino detection as expected in the Pierre Auger Observatory [p|-[TT 



Although CR experiments have detected showers with large zenith angles, most of these data 
have not been yet analyzed due to the difficulties present in understanding the space-time 
structure of the inclined showers. It is worth mentioning that the event reconstruction is 
performed in all air shower arrays assuming that the observation depends only on perpen- 
dicular distance from the shower axis. This very strong assumption, if not satisfied, may 
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introduce important errors in core position and energy determination; precisely this is the 
case for large zenith angle showers. The main difference in analyzing inclined showers lies in 
the fact that the assumed circular symmetry for the observed particle densities from vertical 
showers is broken. Besides, there is an asymmetry due to the different path traveled by 
particles in the upper and lower sides of the plane perpendicular to the shower axis (atten- 
uation effect). These effects which are small at 6' < 45° become important at larger zenith 
angles, producing a significant asymmetry of the lateral distribution of particle density. On 
the other hand, the charged particle trajectories through the atmosphere are distorted in the 
geomagnetic field, effect extensively studied ||l2|,p|,[l3|-p!5|] . A recent analysis of the influence 



of the magnetic field of the Earth on the shower observables for the Southern site for the 
Pierre Auger Observatory (average Geomagnetic Field of 25000nT and site altitude at 1400 
m.a.s.l.) |jl6[ has been performed. The study indicated that deflections alter significantly 
the lateral distribution of muons only at zenith angles > 75°. This result, however, would 
be different if one considers different geomagnetic fields and ground level altitudes ^j. In 
all cases, for smaller angles, the asymmetry in the particle densities is dominated by both 
a purely geometrical effect and the effect on the longitudinal development. As a result, the 
azimuthal symmetry of the spatial distribution is destroyed, changing the dependence on the 
distance from the shower core, directly affecting the estimate of the energy of the primary 
particle. Much of the interpretation of EAS phenomenon involves the LDF; in addition to 
energy determination, the parameter characterizing the actual stage of the EAS develop- 
ment (the so-called age parameter) is expected to contain information about the primary 
composition as well as the characteristic of high energy interaction. 

In this paper, we present an approach to analyze lateral distributions of particles pro- 
duced by UHECR, which is valid at all zenith angles. We tackle the problem studying, by 
means of Monte Carlo simulations, the evolution with atmospheric depth of the parameters 
of semi-analytical lateral structure functions for vertical showers. The results obtained are 
then used to calculate the corresponding LDF for inclined showers including for the first time 
both, geometrical and attenuation effects. We also present a comparison of the proposed 



LDFs with AGASA and Haverah Park experimental data and with MC data including water 
Cerenkov detector effects. 



II. ANALYSIS OF VERTICAL SHOWERS 

A. Electron lateral distribution function 

We start our analysis using the lateral structure function derived for a pure electromag- 
netic cascade by Nishimura and Kamata, and later worked out by Greisen |T^;T^ in the 



well-known NKG formula. This function has been much used in theoretical developments 
(3D extensive air shower simulation) as well as in the description of experimental results 
|1|,|20|. The NKG formula is given by: 



^^""^ rl 27rT{snkg) r(4.5 - 2snkg) \rj V + rj 

21n(-^) 

with Snkg = 3/(1 H and t = patm{z)dz/XQ, where A^e is the total number of 

electrons, eo is the critical energy , Eq is the primary energy, is the Moliere radius and 
Xq is the radiation length in air. The age parameter, Snkg, characterizes the actual stage of 
the EAS development. 

To obtain the lateral distribution of particles it is necessary to consider the multiple 
Coulomb scattering, which determines the characteristic size of the shower front. The lateral 
development in electromagnetic cascades in different materials scales well with the Moliere 
radius |21] given by = ^s^, with the energy scale Es = m.c^(^)°'^. 

It is straightforward to see that the Moliere radius varies with the inverse of the density 
of the medium. 

. Patrajho) _ 9.6gcm-^ 

Tm — Tra\no) -77- — 77-7- l^Zj 

Patm\i>') PatmK,!^} 

In order to study the development of the shower parameters with t, we use Monte Carlo 
simulations of vertical proton induced showers at 10^^ eV generated using aires. This pro- 
gram is a realistic air shower simulation system, which includes electromagnetic algorithms 



2^ and links to different hadronic interactions models. We use in this paper qgsjet ||23 



for nuclear fragmentation and inelastic collision which is consistent with experimental data 
2M. For the highest energy air showers the number of secondaries become so large ( > 10^^) 



that it is prohibitive, in computing time and disk space, to follow and store all of them. 
A. M. Hillas introduced a non uniform statistical sampling mechanism which allows to 
reconstruct the whole EAS from a small representative fraction of secondaries that are fully 
tracked. Statistical weights are assigned to the sampled particles to account for the energy 
of the discarded particles. This technique is called "statistical thinning" . aires includes an 



extended thinning algorithm, which is explained in detail in reference 1^^. The present work 
has been performed using in most cases an effective thinning level eth = Eth/Eprim = 10~^ 
which allows to avoid unwanted fluctuations and provide a reliable number of particles far 
from the shower core. All shower particles with energies above the following threshold were 
tracked: 90 keV for 7's, 90 keV for electrons and positrons, 10 MeV for muons, 60 MeV for 
mesons and 120 MeV for nucleons and nuclei. 

With the particle data recorded we have evaluated the lateral distribution of electrons 
and positrons not only at ground altitude, but also at different atmospheric depths t. Fits to 
lateral distribution functions (LDF) were performed using equation (1) in order to obtain the 
corresponding parameters as a function of t, in a shower-to-shower analysis. We found that 
the Moliere radius deviates from (2) when using a value of the age parameter fixed to the 
theoretical prediction value, which is an approximation valid for pure electromagnetic cas- 
cades. Undoubtedly, a different age parameter should somehow be used to reproduce showers 



initiated by hadrons. This point has already been addressed experimentally p8|- p2| , p4| -|36 
and extensively studied by several authors(See and references therein). 

Our studies of vertical showers confirm, as proposed by |3^, that a modified NKG formula 
should be used. Modifications of the NKG formula have been presented using theoretical 



considerations [^] as well as parameterizations induced by Monte Carlo simulations . In the 
first case, the method is subject to several approximations in order to obtain an analytical 
expression. The second one relies on Monte Carlo results at ground level yielding to average 



values of the parameters which are then included in the parameterizations. For our purpose 
it is important to obtain a very general LDF which can be used to describe the behavior of 
lateral distributions of electrons and positrons for vertical showers, not only at ground level 
but at any atmospheric depth. We use a NKG function with an age parameter given by 

* ^ 1+2/3 /t ("^) 

where (5 takes into account the above mentioned deviations from the theoretical value Snkg- 
The proposed NKG formula provides a good description of the LDF at all stages of shower 
development only outside the hadronic core which lies close to the shower axis. It should 
be mentioned at this point that most UHECR experiments can only measure densities at 
r > 100m from the axis, since in the very dense region near the shower core the particle 
numbers are so large that array detectors would be saturated for high energy showers. Fits 
to Monte Carlo data at r > 100m from the shower core were performed with Vm fixed to (2), 
leaving the overall normalization and /3 as free parameters. 

Figure 1 (top) shows the Monte Carlo e+e~ density distributions corresponding to a 
single 10^^ eV proton shower at selected atmospheric depths from the whole set of levels 
considered. The error bars are included in all cases in the points. Sohd hues are the results 
of the fits. The reduced are in the range 1-3 with better fits around the position of the 
shower maximum. 

The behavior of the fitted age parameter with atmospheric depth is displayed in Figure 
2. A set of 15 proton initiated showers was used to analyze the natural fluctuations. With 
different symbols it is shown the average fitted s value at different distances from the shower 
core; 150 m < r < 450 m (circles), 150 m < r < 750 m (squares) and 150 m < r < 1600 
m (triangles). The error bars, in most cases smaller than the points, represent the natural 
shower-to-shower fiuctuations (The RMS fluctuations of the mean are always smaller than 
the symbols). A remarkable feature is the fact that the whole radial range beyond 150m 
can be described by a unique age parameter independent of r. The band in the same figure 
shows the age parameter calculated using (3) with the assumptions that (3 is given by the 
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position of average shower maximum, < Xmax >= 730 ± 36g cm~^ for the shower set used. 
It is evident that, around the shower maximum the agreement is better and the value of the 
age parameter is approximately 1 at Xmax- 

The total number of electrons obtained from the fit to each single shower is slightly lower 
than the actual one. This can be explained due to the fact that the shower size depends on 
extrapolation of the LDF in the region close to the core, which is not included in the fit. 



B. Muon lateral distribution function 

It is well established that at ground level and for large zenith angles the muonic com- 
ponent of an EAS becomes dominant [ff6|,0, reaching a maximum proportion at 75°, after 



which it decreases slightly |T6|. Consequently, the study of muon density distributions at 
all zenith angles is very important. Several parameterizations have been presented for the 
muon lateral distributions produced in vertical showers; a NKG-type LDF was empirically 



derive by Greisen 



ly \ 0.75 / ■y^ \ 2.5 



f,{r,s)^[-\ 1 + - 



fixing tg = 320m. Vernov et al( and references therein) proposed a semi-analytical form 
of the structure function 

c(i-)-%xp(-^) (4) 

with 7 = 0.4 and tq = 80m. A similar approach was suggested by Hillas The slope of the 
two functions are in very good agreement at intermediate distance, but Vernov distribution 
is flatter close to the shower core and decreases faster at larger distances. These LDF 
have been used to fit experimental data, however, neither function seems to reproduce the 
whole radial range of an EAS. The reason for that is clearly the lack of dependence on 
shower stage development, using instead fixed parameters in the corresponding structure 
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functions. Very recently, the KASCADE experiment has used a NKG formula to fit muon 
density distributions. The fits were performed close to the shower core (r < 200m) with non 



conventional and s values [jl9| 



A dependence of the parameters according to the shower age must be expected for the 
LDF although not precisely the same as the one exhibited by electromagnetic component. 
Unlike electrons, muons in an air shower are less attenuated and little affected by Coulomb 
scattering. Its spread is determined by the direction of emission of the parent particle and 
hence increases while the shower propagates downwards. A shower age dependence of the 



muon structure function has been considered in relation with the parameter 7 (See [0 and 
references therein). 

For the purpose of this paper the Vernov function seems to be a better choice due to 
its analytical origin. In order to study the depth dependence of the parameters, we use a 
modified Vernov distribution with a 7 parameter dependence with atmospheric depth given 
by 7 = 2 — s, with s as in (3) due to the fact that in the radial region smaller than the scale 
radius, the slopes predicted by both the NKG-like and Vernov distribution are very similar. 

It has been shown pT| , ^ that, if the parent particles are created with a p± distribution: 
P±/Po^wi~P±/Po)dp±/po, then the Vernov distribution is obtained at ground level, with tq 
given by 

where < p± >= 2pQ is the mean transverse momentum, < E^> the mean energy of muons 
and < Hp > the mean height of production. 

These approximate expressions can serve to calculate the variation with t of the param- 
eters characterizing the lateral spread. The ratio < p'^ > / < > can be considered 
constant while the shower develops [0 and the variation of tq with altitude is only deter- 
mined by the t dependence of the mean high of origin of muons. These particles are produced 
in every pion generation and their energy distribution follows that of their parents. Thus, 
one expects the dependence of < Hp > with t to arise from two sources. The trivial one 



resulting from the fact that we consider observing levels at different altitudes, which gives a 
factor ho \og(t/tp) assuming an isothermal atmosphere of scale height ho; and a dynamical 
one due to the pion decay process. It is well known that, at a given energy, the pion decay 
probability decreases with increasing air density. Thus, a characteristic radius ro{t) results, 

where to is the ground level depth in gcm"^. Note that aX t — to equation (5) is 
recovered. 

Fits using the modified Vernov function to the muon density distributions obtained from 
the same Monte Carlo set of showers generated as described in the previous section, were 
performed with A^^, ro, and /3 simultaneously free. Reduced of the fits to single showers 
for all predetermined observing levels are in the range 1-3. Figure 1 (bottom) displays the 
^JL~ density distributions (points) for a single 10^^ eV proton shower corresponding to a 
subset of the atmospheric levels, together with the fits (sohd hues) where the agreement is 
evident. 

The parameter (5 is approximately constant and much bigger than the shower X^ax. 
This feature is somehow consistent with the fact that muons become relatively more and 
more abundant as the altitude decreases with only a slight attenuation. 

The mean value of Tq from the above mentioned set of showers is plotted in Figure 3, 
where error bars indicate the natural shower-to-shower fluctuations. The solid line corre- 
sponds to a fit using equation (6) which gives tp of the order of A^v, the attenuation length 
for nucleons in air, as expected. Note that ro is a flat function of t resulting from a com- 
pensation of the two factors described above. For inchned showers, the typical slant depth 
of the shower increases significantly with the zenith angle and the t dependence becomes 
dominated by the atmospheric density factor. 

Figure 4 displays the obtained age parameters fitting different radial ranges: 100-300 m 
(circles), 100-600m (squares) and 100-1600 m (triangles). The error bars indicate as in the 
other cases the natural shower fluctuations. The curve shows the prediction using (3) and 
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constant value of /3 = 1800 g cm~^. The total number of muons A*"^ from the fits agrees 
quite well with the corresponding values registered from the Monte Carlo data, in spite of 
the fact that the fits are performed in all cases at core distances r > 100m, implying an 
extrapolation of the LDF beyond the range considered. 



III. INCLINED SHOWERS 

The main goal of the analysis done on vertical shower was to obtain the atmospheric 
depth dependence of the shower parameters. This result is crucial to properly take into 
account the effect of the attenuation in the atmosphere observed in inclined showers, as we 
shall show in what follows. In this section we extend the solutions obtained for vertical 
showers to inclined ones, considering not only the geometrical effect but also the influence 
of the attenuation in the atmosphere. 

It must be emphasized that the measured data by surface arrays are ground density of 
particles and the showers are analysed projecting the ground densities into the shower plane. 
As mentioned in the introduction, this procedure allows correction of the geometrical effects 
and it is a good approximation for near-vertical showers. However, it was shown in the 
section II for vertical showers, that the parameters of the LDF depend on the atmospheric 
depth t. This means that to explain the observed asymmetry, besides the transformation to 
shower front coordinates, one has to evaluate all parameters of the LDF at the new slant 
depth t' 

^, ^ !T Patm{z)dz' 

with z' along the shower axis direction, in order to include the attenuation effect in the 
atmosphere. 



Actually, the structure of an air shower is quite complex but following |P3| , |4^ , an inverted 
cone seems to be a reasonable first approximation. The intersection of an inclined cone with 
a plane parallel to the ground gives. 
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— L — (z' — L) 1 tan 6 cos (f) 

cos 9 \ z— L 

The parameters L, 6 and are indicated in Figure 5. If the cone angle a is assumed 

nearly independent of z then, the slant depth t' results from: 

— = cos^^ (1 — tana tan^^ cos^') (9) 

t 

Particle density distributions for inclined showers are given by the corresponding electron 
and muon LDF but evaluated at slant depth t' where the dependence on the azimuthal angle 
is evident. Proton air showers at 10^^ eV were generated under the same conditions of vertical 
showers with AIRES + QGSJET at zenith angles between 9 — 10° and 9 — 75°, where the 
influence of the geomagnetic field is still negligible. The geomagnetic field was then switched 
off in the simulation. The primary particle arrives at the corresponding zenith angle from 
0° East of the South. 

For single showers, two-dimension density distributions r-0 were built at all predeter- 
mined observing levels for different particles with bins 100 m and 30 deg wide. 

A. Electron lateral and azimuthal distributions 

The distributions were fitted using the modified NKG formula with s (Eq.(3)) including 
the complete zenith angle dependence, with a, (5 and as free parameters. The proposed 
LDF fits very well e+e~ distributions at all observing levels for 9 < 60°. At zenith angles 
greater than 60° the fits arc good only up to a slant depth t' ~ 2000g cm~^. This can be 
understood at the light of the physics behind air showers induced by nucleons or nucleus. 
These particles, undergo nuclear coUisions in the atmosphere. The cascade begins with 
a hadronic interaction with an increasing number of pions in each generation of particle 
interactions. The rapid decay of 7r°'s into photons feeds the electromagnetic cascade, which 
finally dissipates ~ 90% of the primary energy while charged pion decays produce the muons 
and neutrinos. The development of the shower reaches a maximum at around 800 g cm~^ 
for vertical 10^^ eV proton showers. However for inclined showers, the slant depth rises 
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dramatically from ^ 1000 {6 = 0°) to 36000 g cm'^ {6 = 90°) . A detailed study of the 
characteristics of the muon and electromagnetic component in very inclined showers can be 
found in reference [0]. A standard feature of very inclined showers is that particles at ground 
are dominated by the muonic component, while the electromagnetic component is almost 
completely absorbed due to the greatly enhanced atmospheric slant depth. Beyond ~ 2000 
g cm~^ the small electromagnetic cascade is mostly originated from muon decay and one 
cannot expect a NKG distribution for electrons originated from muon decay. 

Figure 6 shows a set of r projection (integrated contents over all 0) and projections 
( integrated contents over the fitted r range) on the shower front from the two-dimension 
distributions at given atmospheric depths, corresponding to 30° and 60° inclined showers. 
Points with error bars stand for simulated data and curves represent the result of the fit. It 
is important to note the asymmetry due to atmospheric attenuation effects, higher particle 
density at 0° azimuth angle and a depression at 180°. This is evident in the plots showing 
(f) projections, even for the case 6 = 30°, becoming very important at larger zenith angles. 
Note also the symmetry in the density at ±90°. It can be seen from the plots that our 
approach for the LDF predicts the correct dependence for the azimuthal angle distribution. 
In observations with large zenith angles, the position of the depth of maximum is more 
distant from the ground with an increasing atmospheric absorption of the electromagnetic 
component. However, the plots indicate that, even at 6 < 70° a significant number of 
electrons can reach the ground level. 

Figure 7 (top) displays the fitted age parameters at selected atmospheric depths for 
single showers at 0° (circle), 30° (square), 50° (triangle) and 60° (diamond). When plotting 
together the results obtained from different showers at different levels, it is put into evidence 
that the fitted age parameters overlap smoothly with only small differences due to shower- 
to-shower fluctuations. This result implies a consistency check and confirms the validity of 
our approach. With the same aim, the corresponding total number of electrons obtained 
from the AIRES output are plotted in solid lines in Figure 7 (bottom), along with the fitted 
shower size Ne (points with same convention than in Figure 7 (top)) at selected atmospheric 
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depth. It is once more worth mentioning that the curves which correspond to single showers 
with different incident zenith angles are of course, affected by the fluctuation of the first 
interaction point. However, the shower size obtained from the fit agrees rather well with the 
values from the Monte Carlo results, showing the expected longitudinal shower development. 

The parameter tan(Q;) takes values between 0.15 and 0.2, as expected, since an EAS 
covers a large area but subtends a small apex angle. 

B. Muon lateral and azimuthal distributions 

The same set of generated showers, proton air showers at 10^^ eV at zenith angles between 
9 = 30° and 6 = 70° with an effective thinning level of 10~^, was used to analyze the muon 
component. 

Using the procedure described in the previous section for electrons, for each single shower 
two-dimension density distributions r-0 were built at all predetermined observing levels. In 
order to perform the fit to the muon densities resulting from the inclined showers, the mod- 
ified Vernov LDF proposed in section II. B was used. The azimuthal angle dependence was 
included following the approach described in section III. The proposed function successfully 
fits the Monte Carlo data at all zenith angles. Reduced are in the range f^i 1 — 3. As in the 
case of electrons, the fitted A*"^ and s values are compatible with the expected ones. Figures 
8 show the radial and azimuthal angle projections of the muon density from two-dimensional 
fits to the data for 60° and 70° zenith angles. The dots stand for the Monte Carlo data and 
the curves correspond to the fits. It seems worth noting that the asymmetry introduced by 
the atmospheric attenuation in the lateral distribution is small, mainly at observing levels 
near the ground. This feature contrasts with the observed one for electrons, even at lower 
zenith angles. In previous analysis of quasi- horizontal showers 0, the attenuation effect was 
neglected. However, as mentioned in the introduction, for showers arriving at > 75° (with 
high /i content at ground level) the lateral asymmetry is dominated by geomagnetic effects. 
The complex behavior of the muon component due to magnetic field effects was studied in 
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great detail in reference [7] . 



IV. COMPARISON WITH EXPERIMENTAL LDF 

As mentioned above, to estimate the energy of the primary particle, analysis of the lat- 
eral distribution of shower particles at ground are important. Depending on the particle 
detectors of ground array experiments, different LDF parametrizations have been used. The 
signal measured by plastic scintillators (Volcano Ranch, Yakutsk and AGASA) is dominated 
by the electromagnetic component of the shower compared with the muon signal. In Water 
Cerenkov detectors (as used in Haverah Park and being implemented for the Pierre Auger 
Observatory) the electromagnetic and muon components contribute to the total signal. Fig- 
ure 9 shows the signal produced by different particles (top) and the ratio of the muon to 
electromagnetic Cerenkov signal (bottom) measured in units of vertical equivalent muons 
per m^ for a 10^^ eV vertical proton shower. It is evident from the plot that the signal 
changes with distance from the core, being dominated by the muon component far from the 
shower axis. 

The implementation of the method proposed in this paper to obtain the LDF for inclined 
showers once the appropriate parametrization for vertical showers is found, clearly implies 
the analysis of the parameters of the corresponding LDF to find their dependence with 
depth. This specific study has to be performed in all cases using simulations. 

In what follows, we compare the parametrizations proposed in section II and III for 
particle densities with experimental data. It should be stressed that we do not intent to 
re-analyse the experimental data but rather to investigate to which extent the obtained LDF 
from pure shower simulations are valid to reproduce experimental data without including 
additional modifications. In the case that the LDF are suitable to describe experimental 
data, it would be possible to further check the method presented in section III. 

For scintillator detectors, generahzations of the NKG distribution have been used to 
parametrize the observed density with distance. In Figure 10 we present the experimental 
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result of AGASA for vertical 10^^ eV showers [B5| in the form of their empirical functions 



for the total signal S(r) and density of muons of threshold energy 1 GeV, together with 
the result of the fit using the parametrizations proposed in this paper. The S(r) signal was 
fitted with the modified NKG function of section II, with the conventional assumption of 
Tm = 91.6 m. This function can reproduce the AGASA experimental LDF over the core 
distance 100-2500 m surprisingly well with a unique value of the age parameter s = 1.05. If 
one assumes a pure electromagnetic cascade, it is possible to make a very crude estimation 
of the primary energy from the fitted f3 parameter which results to be E' = 1.7 x lO^^eV. The 
value of the signal density at 600 m from the core with our fit is 41 m^^ to be compare with 
the one obtained from the AGASA results S'(600) =51 m~^. The average conversion relation 
for primary proton and iron using two different hadronic interaction models determined by 
the AGASA experiment is £"[6^] = 2.23 x 10^'^ 5*0(600)^ °^ [0. Using this expression we can 
estimate the energy from our fit to he E = 1.2 x lO^^eV. 

In addition, we fitted the AGASA muon lateral distribution with the modified Vernov 
function from section II. B in the range 100-2500 m with all parameters left free. The result 
of the fit yields to tq = 436 m. It can be seen that the shape of the fitted LDF is consistent 
with the empirical AGASA function. The value of p^(600) from our parameterization is 3.86 
m~^ in perfect agreement with the value obtained using the AGASA function. 

Due to the lack of published signal densities in the X-Y plane for large zenith angle 
showers from AGASA experiment, we present here only comparisons with vertical showers. 
It would be straightforward to extend the analysis to non-vertical showers using the same 
LDF and the method we proposed in section III. It is worthwhile to remark that most of 
the experiments have only published cosmic ray showers at small zenith angles. 

Water Cerenkov signal (WCS) due to EAS particles at ground includes both electro- 
magnetic particles and muons, and the signal produced in the detectors is the sum of both 
components. In Figure 11 we show an extremely large EAS observed by the Haverah Park 
experiment and their predicted LDF in the regions r < 800m and r > 800m as described in 
(solid lines). The result of the fit using the modified NKG formula from section II. A is 
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also plotted. The Moliere radius was fixed to the value calculated using (2) and the mean 
atmospheric depth of the array; 1018 g cm~^. Our LDF reproduces rather well the Haverah 
Park data in the whole range with age parameter s = 1.24. From the fitted function a 
value of p(600) = 23.8 VEM m~^ was obtained to be compared with the value obtained by 
Haverah Park p(600) = 25.2 VEM m'^. 

In order to check the validity of our approach for Water Cerenkov detectors we have also 
made estimations of the WCS using a direct conversion procedure that retrieves average 



signals for every particle hitting the detector [^. Such averages were evaluated using the 
AGASIM program ||4^. Fig. 12 (Top) shows the obtained WCS for all the particles at 
ground level for a 10^^ eV vertical shower for distances to the axis beyond 200 m (dots). 
The muon threshold energy is 300 MeV. The result of the fit using the modified NKG formula 
leaving all parameters free is superimposed as solid lines, providing a good description of the 
data {x'i/d.o.f = 26/12 ^ 2). The scale radius determined by the fit is 1.4 times the pure 
theoretical Moliere radius. The age parameter results s = 1.15, being strongly correlated 
with tq. From the fitted function a value of p(600) = 12.7 VEM m~^ was obtained. With 
this value and the energy calibration formula used to analyse the Haverah Park data ED 
{E{EeV) = 0.612 x p(600)°-^^), the estimated primary energy results Ep ^ 0.8 x lO^^eV. 
Using the parametrization for non-vertical showers derived by application of our method to 
include the effect of the attenuation in the atmosphere in the modified NKG function, we 
also fitted the lateral density of WCS at ground for a proton shower at zenith angle of 40°. 
At this modest zenith angles the average signals estimated as described above are still the 
appropriate ones. The resulting LDF fits well the data {x^/d.o.f = 383/165 = 2.3) with a 
fixed value of tq = lAvm as for the vertical showers. The s parameter tends to decrease. 
It is, however, well known that the shower parameters s and tq are highly correlated. The 
fitted parameter tan(a;) takes a value of 0.1 in agreement with the expected one. In Figure 
12 (bottom) the radial and azimuthal WCS density projections, measured in VEM m~^, are 
plotted together with the result of the fitted LDF. 
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V. CONCLUSIONS 



It is clear that information concerning the lateral distribution functions of particles is 
very important to determine the characteristics of the ultra high energy cosmic rays. In this 
paper, we presented a new approach to obtain appropriate parameterizations of LDF, able 
to reproduce particle density distributions generated for UHECR arriving from any zenith 
angles. 

We started our study with the analysis of vertical showers using modified NKG (section 
II.A) and Vernov (section II. B) structure functions for electrons and muons, respectively. 
A detailed study of the shower parameters with atmospheric depths was performed. The 
results indicate that a unique age parameter can be used at radial distances greater than 150 
m from the shower core, which is the region of interest for surface arrays. The knowledge 
of the behavior of shower parameters as the cascade of particles develops is crucial to work 
out the LDF for inclined showers. It is worth noticing that the method presented here 
it is well suitable to a shower-to-shower based analysis. We paid special attention to the 
fact that the LDF of particles from inclined cascades must be corrected not only due to 
geometrical effects but also for the longitudinal attenuation. This latest, usually neglected 
or considered in a very simple manner in the experimental analysis, is surprisingly large for 
electron distributions, even at moderate zenith angles where the amount of electrons reaching 
ground level dominates. Thus, data analysis assuming that the lateral distributions depend 
only on the distance from the shower axis might yield to biased results if the atmospheric 
attenuation effect is not considered at all. Muon distributions are also affected by attenuation 
effects, although to a less extent. The asymmetry at very large zenith angles in this case is 
dominated by geomagnetic effects. 

The obtained lateral distribution functions fit very well the generated data at all observ- 
ing levels along the longitudinal shower paths. The very good agreement of the theoretical 
distributions for inclined showers obtained using our method with the radial, and more 
important, the azimuthal angle distributions, provides strong supporting evidence for the 
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validity of this approach. In connection with this result we have checked our method with 
the signal detected by ground array detectors. We have shown that if the signal produced in 
water Cerenkov tanks is parametrized by a modified NKG function, then the application of 
our method is straightforward and the corresponding LDF for inclined showers reproduces 
the experimental density distributions rather well. 

A final comment on the method proposed is well worth making here: Several present 
experiments use both NKG and Vernov type LDF. The application of our method would 
allow them to include high zenith angle showers in their analysis, what implies a significant 
increase of their acceptance. Future experiments will also benefit from the analysis of showers 
arriving from any direction 
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FIGURES 




FIG. 1. 



Fig. 1: Electron (top) and muon (bottom) lateral distributions of a 10^^ eV vertical 
proton shower at different atmospheric altitudes. The solid line indicates the result of the 
fit using the proposed LDF (See text). 
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FIG. 2. 

Fig. 2: Atmospheric depth dependence of the average age parameter obtained fitting 
electron lateral distributions of single showers at different radial distances: 150 m < r < 450 
m (circles), 150 m < r < 750 m (squares) and 150 m < r < 1600 m (triangles). Shaded area 
indicates the expected dependence using a NKG formula with (1). 
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FIG. 3. 



Fig. 3: Atmospheric depth dependence of the average scale radius ro(t) fitting single 
showers. The error bars correspond to the standard fiuctuations. Solid line is the result of 
the fit using equation (6). 
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FIG. 4. 

Fig. 4: Atmospheric depth dependence of the average age parameter obtained fitting 
muon lateral distributions of single showers at different radial distances: 100 m <r < 300 
m (circles), 100 m < r < 600 m (squares) and 100 m < r < 1600 m (triangles) The line 
indicates the fit as described in the text. 
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FIG. 5. 



5: The intersection of an inclined cone and the ground plane. 
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FIG. 6. 



Fig. 6: Projections of two-dimension electron distributions r-0 at different slant depths 
for a lO^^eV proton shower at 30° and 60° zenith angles. 
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FIG. 7. 



Fig. 7: Upper plot: Age parameters at selected atmospheric depth for showers entering at 
0° (circle) ,30° (square), 50° (triangle) and 60° (diamond) zenith angles. Shower size obtained 
from the fits vs slant depth are shown in the lower plot. The actual total particle numbers 
for each single shower are indicated with solid lines. 
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FIG. 8. 



Fig. 8: Projections of two- dimension muon distributions r-(j) at different slant depths for 
a lO^^eV proton shower at 60° and 70° zenith angles. 
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FIG. 9. 



Fig. 9: Top: Electron (open circles), gamma (squares), muon(triangle) and total (circles) 
water Cerenkov signal in VEM m~^ vs distance from the shower core for a 10^^ proton vertical 
shower. Bottom: Ratio of muon to electromagnetic signal vs distance from the core. 
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FIG. 10. 



Fig. 10: Lateral distributions of total signal (solid line) and muon(dotted line) deter- 
mined by AGASA for vertical lO^^eV showers, compared with fits using modified NKG 

(dashed line) and Vernov (dash-dotted line) LDF. 
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FIG. 11. 



Fig. 11: EAS observed by Haverah Park. The solid line indicates the Haverah Park 
predicted LDF and dashed-line corresponds to the result of the fit using the modified NKG 
LDF. 
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FIG. 12. 



Fig. 12: Top: Lateral distribution of signal produced by electrons, muons and photons 
in a water Cerenkov detector in units of vertical equivalent muons (see text). Bottom: 
Projections of two-dimension water Cerenkov signal distributions r-0 for a proton 40°. Solid 
line corresponds to the fit performed usig the modified NKG LDF. 
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